Stationary entanglement induced by dissipation 
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The dynamics of two two-level dipole-dipole interacting atoms coupled to a common electromag- 
netic bath and closely located inside a lossy cavity, is reported. Initially injecting only one excitation 
in the two atoms-cavity system, loss mechanisms asymptotically drive the matter sample toward 
a stationary maximally entangled state. The role played by the closeness of the two atoms with 
respect to such a cooperative behaviour is carefully discussed. Stationary radiation trapping effects 
are found and transparently interpreted. 

PACS numbers; 03.65.Yz, 03.67.Mn, 42.50.Fx 



A system prepared in an entangled state may show 
rather puzzling behavior and counterintuitive features. 
In such a state, for example, the system exhibits the as- 
tonishing property that the results of measurements on 
one subsystem cannot be specified independently from 
the parameters of the measurements on the other com- 
ponents. The renewed and more and more growing in- 
terest toward entanglement concept reflects the consol- 
idated belief that unfactorizable states of multipartite 
system provide an unreplaceable applicative resource, 
for example, in the quantum computing research area 
0. However, the realization of quantum computation 
protocols suffers from the difficulty of isolating a quan- 
tum mechanical system from its environment. Unavoid- 
ably, the system-environment interaction leads to deco- 
herence phenomena which, as intuition suggests, are al- 
ways noxious for quantum computers, since they imply 
loss of the information stored in the system. This cir- 
cumstance is at the origin of an intense research aimed at 
proposing theoretical and experimental schemes to fight 
or control decoherence manifestations IE ^ IE IE IE III ■ 
Very recently, however, the possibility of environment- 
induced entanglement in an open quantuni system has 
opened new intriguing perspectives [E 0, ^| ■ For 
example, it has been shown that dissipation can be ex- 
ploited to implement nearly decoherence-free quantum 
gates ^E J the main requirement being the existence 
of a decoherence-free subspace for the system under con- 
sideration 0. Moreover, despite the widely held be- 
lief, it has been shown that transient entanglement be- 
tween distant atoms can be induced by atomic sponta- 
neous decay or by cavity losses |0, IT^ 0^ ■ It has been 
also demonstrated |[17| that asymptotic (stationary) en- 
tangled states of two closely separated two-level atoms 
in free space can be created by spontaneous emission. 
Notwithstanding the nearness of the two atoms, dipole- 
dipole interaction is however neglected in Ref. [TE|- In 
this paper we investigate the dynamics of a couple of 
spontaneously emitting two-level atoms confined within 
a bad single-mode cavity, taking into account from the 
very beginning their dipole-dipole interaction. Our main 



result is that in such a condition the matter subsystem, 
even experiencing a further relatively faster energy loss 
mechanism, may be as well conditionally guided toward 
a stationary robust entangled state. We shall prove that 
this state is the antisymmetric one with respect to the 
exchange of the two two-level atoms provided that all 
the initial energy given to the system is concentrated in 
the matter subsystem only. 

Let us suppose that two identical two-level atoms are 
located at ri and r2 inside a resonant single mode cavity. 
Putting R = |ri — = |R| and denoting by 6 the an- 
gle between R and the atomic transition dipole moment 
d, the Hamiltonian describing the dipole interaction be- 
tween the two atoms can be written in the form ,18> .1^ 

Hi2^hT^{a^_^^<j'--^ +h.c.), (1) 

where rj = | ,^°j^3 (1 — -icos^O) and Fq is the spontaneous 

emission rate in free space. In eq. {i = 1,2) are 

the Pauli operators of the i-th atom. Assume in addition 
that all the conditions are satisfied under which the inter- 
action between each atom and the cavity field is described 
by a Jaynes Cummings (JC) model 20]. Under these hy- 
potheses, the unitary time evolution of the system under 
scrutiny is governed by the following Hamiltonian 

.2 2 

Hac = Wa-f^^a«+;i^[e«aa^'^-t-/i.c.]+Fi2, 

i=l i=l 

(2) 

where Iuvq denotes the energy separation between the 
i-th atom [i = 1,2) excited (|+)i) and ground 
states, uj ujQ \s the frequency of the cavity mode 
and e*^*' is the i— th atom-field mode coupling constant. 
In the above equation a {a^) is the annihilation (cre- 
ation) operator relative to the single-mode cavity field 
and iJz the i— th atom inversion operator. It is easy to 
demonstrate that the total excitation number operator 

N = a^a -\- \{o''z'^ + o-f^) -f 1 is a constant of motion. 
Thus, preparing the physical system at i = in a state 
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with a well defined number of excitations N^, its dynam- 
ics is confined in the finite-dimensional Hilbert subspace 
singled out by this eigenvalue of N . In a realistic situa- 
tion, however, the system we are considering evolves un- 
der the action of different sources of decoherence. First, 
photons can leak out through the cavity mirrors due to 
the coupling of the resonator mode to the free radiation 
field. Moreover, the atoms confined in the resonator can 
spontaneously emit photons into non-cavity field modes. 
The microscopic Hamiltonian which takes into account 
all these loss mechanisms can be written in the form l21l 



and 



H — Hac + Hr + Har 
where Hac is given by Eq. (jSJ, 

Hr 



H 



CR, 



is the Hamiltonian of the environment, 



HAR^J2^9^Z£^,x4^+h.. 



(3) 



(4) 



(5) 



describes the interaction between the atoms and the bath 
and, finally. 



HcR = ^[sk,ACk,Aa^ + h.c] 



(6) 



k.A 



takes into account the coupling between the environment 
and the cavity field. As usual, in Eq. we have assumed 
that the two subsystems, the two atoms and the single- 
mode cavity, see two different reservoirs in the following 
both assumed at T = 0. The boson operators relative 
to the atomic bath are denoted by {ck,A, cj^ x} whereas 

cu,\,cl^x are the (k. A) mode annihilation and creation 
operators respectively of the cavity environment. More- 
over, the coupling constants Sk,A are phenomenological 
parameters and 



ffk,A 



-li -)-^ekx ■ 



(7) 



stems from a dipole atom field coupling In Eq. (0) 
GkA represents the polarization vector of the thermal bath 
(kA) mode of frequency Wk and V its effective volume. 
In the Rotating- Wave and Born-Markov approximations 
|2^ |23 | , the reduced density operator pAc relative to the 
bipartite system composed by the two atoms subsystem 
and the single-mode cavity, evolves nonunitarily in ac- 
cordance with the following quantum master equation in 
Lindblad form: 



PAC = -j^[Hac,Pac] + C,fPAc + Cap AC, (8) 



where 



C-fPAC = k{2apAca^ - PAca - pAca^a), (9) 



CaPac — 



(10) 



^T{2(j\' pAC^X - (y\ c!- PAC - PAco\'a\') + 

i=l 

-p /o (1) (2) (1) (2) (1) (2)n , 

Vx2\20- PAC^X ~ oX'a^' PAC - PAC(^+ ) + 

-P /o (2) (1) (2) (1) (2) (1)n 

T 2l[2(j\ ' PAC - (T\'(7'_' PAC - PAC(^+ )■ 

In the above equations we have introduced T = ^^^^jf 
and = T21 = Tf{R) where 

fiR) 



3,. . .9, sinf^i?) 

Z OJo-tl 



(11) 



,,,, „cos(^R] ^sm(^R). 



}■ 



with < f{R) < 1. We emphasize that mea- 
sures the strength of the atom-atom cooperation induced 
by their coupling with a common bath. Finally, k = 
EkA l'SkAp<5(wk — w), appearing in Eq. lO, is the cavity 
decay rate coefficient. Assume now that the two atoms 
inside the cavity are closely located that is c << Rojq. In 
such a situation the cooperation between the two atoms 
stemming from their interaction with a common bath, is 
maximum that is f{R) = 1. Moreover it is reasonable to 



put e 



(1) 



r(2) 



£. We are able to exactly solve Eq. 



when this point-like model is adopted. 

To this end consider the unitary operator U defined as 



C/ = exphJ(a>^-aia^)]. 



(12) 



It can be shown that [U, N] ~ and that, confining our 
attention on the Hilbert subspace correspondent to = 
0, 1, Eq. JHJ can be equivalently cast in the form 



PAC = U'' pAcU ^ --[Hac, PAc] 
n 

+ k{2apAca^ - a^apAc - PAcct^ct) 



(13) 



OT-^O (1)- (1) (1) (1) ~ ~ (1) (1)^ 

2T{2a'_' pac^+ - (^+ Pac - Pac(^+ ) 



where pAc = U^pAcU 

2 

Hac =U^HacU ^ hcoa^'a + hujo^cri'^ + (14) 

1=1 

fee//K^^+/^.c.]-^(af -aW), 

with Eeff = V2e. It is important to notice that in 
the new representation the dipole-dipole interaction H12 
given by Eq. J^l, only renormalizes the atomic frequen- 
cies. The transformed Hamiltonian Hac puts transpar- 
ently into evidence that the system of two atoms cooper- 
ates in the interaction with the cavity field. In this new 
representation two two-level fictitious atoms appear, only 
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one of them being coupled, by means of a simple JC in- 
teraction model, to the cavity field. The circumstance 
that the atomic sample can exchange energy with the 
field through only one of its collective atoms, provides a 
natural explanation for the radiation trapping phenom- 
ena [25l |. This conclusion is true under ideal conditions 
when neither cavity losses or atomic spontaneous decay 
are considered. Thus, it seems interesting to investigate 
whether such an energy storage mechanism survives in a 
more realistic situation like the one under scrutiny in this 
paper. Looking at Eq. H14|l we may see also that in the 
transformed representation the atomic subsystem looses 
its energy only through the interaction of the first atom 
with both the cavity mode and the environment. Such a 
behaviour stems from the fact that the other atom freely 
evolves being decoupled from either the cavity field and 
the electromagnetic modes of the thermal bath. Compar- 
ing this result with the interpretation of the energy trap- 
ping given in Refs. (25i. >26ii] . we may immediately catch 
the main role played by the closeness of the two atoms 
in our model. It is indeed just this feature which, in the 
transformed representation, leads to the existence of one 
atom immune from spontaneous emission losses and, at 
the same time, decoupled from the cavity mode. Thus, to 
locate the atomic sample within a linear dimension much 
shorter than the wavelength of the cavity mode, intro- 
duces a permutational atomic symmetry which is at the 
origin of a collective behaviour of the two atoms. As a di- 
rect consequence the matter system may stationarily trap 
the initial energy even in presence of both the proposed 
dissipation channels. Suppose that only one excitation 
is initially present in the atomic subsystem whereas the 
cavity is prepared in its vacuum state. From an exper- 
imental point of view it seems reasonable to think that 
the excitation given to the matter system is captured by 
the atom 1 or by the atom 2 with the same probability. 
On the other hand, exciting only an a priori fixed atom or 
preparing a quantum coherent superposition of the two 
states |+)i|— )2 and |— )i|-l-)2, could present serious prac- 
tical difficulties. Stated another way, our initial condition 
can be reasonably expressed by pAci^) = Pa{0) (X" pc(0) 
with PA = i(|+)ih)2i(+|2(-| + |-)i|+)2i(-|2(+|) and 
Pc{0) = |0)(0|, |0) being the vacuum Fock state of the 
cavity mode. It is evident that, under this hypothesis, at 
a generic time instant t, the density operator pAC, can 
have non zero matrix elements only in the Hilbert sub- 
space generated by the following ordered set of four state 
vectors 



{|o>|- 



II) 



-)2}. 



(15) 

Taking into account that [U, N] — 0, the same conclusion 
holds for pAc too. This fact provides the key to solve 
exactly eq. (|13ll . One finds 



PAc{t) 



(Pii{t) 

P2.2{t) p2A{t) 

P3,3W 

V p*^^{t) p4.4(i)y 



(16) 



where P3,3(i) = hA^)^ Pii{t) = 1 - Y.i=2 Pai^) '^ith 

+(^^ - lAp) cos(f]2t) + ^) 

sin(02i) + (^ry + ^A_) sinh(17ii))} (17) 



.2 



(cosh(r2i<) — cos{rt2t))e 



-A+t 



(18) 



P^^'^^^ = ffTTI^^^"' .f^i)(sin(f^2i) + 
-zsinh(f7ii)) + Aee//(cosh(f2ii) - cos(f22t))} (19) 



where A±^k± 2r, 2A = -|- iA^ 



(20) 



with P = + 'islff — A"^ and V = —A'^rj'^ and uj — ujq- 
We emphasize that, on the basis of the block diagonal 
form exhibited by pAc{t)i the transformed system, at a 
generic time instant t, is in a statistical mixture of the 
vacuum state |0)|— )i|— )2 and of a one excitation appro- 
priate density matrix describing with certainty the stor- 
age of the initial energy. Taking into account the easily 
demonstrable inequality Oi < it is immediate to con- 
vince oneself that for t ^ A']^^ the correspondent asymp- 
totic form assumed by pAC is time independent and such 
that the probability of finding energy in the effective JC 
subsystem exactly vanishes. Considering that the ini- 
tial condition imposed on pAci^) may be converted into 
Pit = ^{5i2 + Sis) and pij = for i 7^ j, we may conclude 
that at i ^ A'^^ the only two matrix elements different 
from zero are pii{t) = p22{0) and pssit) = ^33(0). Trans- 
forming pAc(t) back to the original representation, the 
exact solution pAC for the reduced density matrix of the 
system under scrutiny is then easily found. Since the uni- 
tary operator U is time independent, we are legitimated 
to forecast an asymptotic time independent solution in 
the original representation too. The reduced density ma- 
trix can indeed be written in the compact form 



PAC = ^|V'r>)(?AD| + ^\i^T){ipT\ 



with IV^d) = |0)|-)i|-)2 and 



(21) 



(22) 



eigenstate of (cti + (72)^ with eigenvalue zero. Eq. l(?T)) 
suggests that stationary entangled states of the two 
atoms can be generated putting outside of the cavity 
single photon detectors allowing a continuous monitor- 
ing of the decay of the system through the two possible 
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FIG. 1: Conditional Concurrence C in correspondence to 77 = 0.5e, 
k = and F = IQ-^e 

channels (atomic and cavity dissipation). Reading out 
the detectors states at f 3> w k^^, if no photon 

has been emitted, then as a consequence of this measure- 



ment outcome, our system is projected into the state \4't) 
given by Eq. (|22|) . This is the main result of our paper 
which means that a successful measurement, performed 
at large enough time instants, generates an uncorrelated 
state of the two subsystems atoms and cavity, leaving the 
atomic sample in its maximally antisymmetric entangled 
state (|22|l . To analyze the time evolution of the degree 
of entanglement that gets established between the two 
initially uncorrelated atoms, we exploit the concept of 
concurrence C first introduced by Wootters [23, ■ If 
at an assigned time t no photons have been emitted, the 
conditional concurrence C assumes the form 

C{t) = I^^^W-^33(t)| (23^ 

P22(t) + P33(f) + pAA{t) 

As clearly shown in Fig. ^ the degree of entanglement 
between the two atoms, starting from zero, reaches its 
maximum value, that is C = 1, in a time of the order of 
k~^, after a large number of oscillations. 

Summarizing, in this paper we have analyzed the dy- 
namics of a realistic system composed of two two-level 
atoms with dipole-dipole interaction, embedded in a bad 
single mode cavity and coupled to a common electro- 
magnetic environment. The exact analytic solution of 
the Markovian dissipative dynamics reveals that the en- 
vironment induces stationary entangled states of the two 
atoms starting from a realistic preparation of the matter 
subsystem. This fact makes it possible the experimental 
realization of this process with real atoms and could be of 
some relevance for the development of quantum devices. 
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